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ABSTRACT 


A  simplified  model  1b  proponed  for  the  combustion  of  a  bi¬ 
propellant  spray.  The  model  considers  two  subsystems:  one  the  spray, 
the  other,  combusted  gases.  These  subsystems  are  coupled  together  by 
heat  transfer,  mass  transfer  and  momentum  transfer.  The  solution  of 
the  problem  relieB  upon  the  integration  of  a  system  of  Beven  nonlinear 
differential  equations.  The  integration  is  performed  using  the  IBM  JOk 
digital  computer.  It  appears  that  under  the  given  boundary  conditions 
a  co-operative  evaporation  process  occurs,  in  which  increased  velocity 
gradients  cause  increased  evaporation  which  increases  the  velocity  gradi¬ 
ent. 


A  series  of  tests  have  been  made  on  a  two  inch  diameter,  variable 
length  rocket  motor,  using  JP-5A  and  liquid  oxygen,  f  simple  converging 
diverging  nozzle  is  employed  to  give  high  chamber  exit  Mach  numbers 
(.4  to  .5).  For  a  given  set  of  injection  parameters,  under  stable  oper¬ 
ating  conditions,  the  results  indicate  that  the  steepest  portion  of  the 
gradient,  which  corresponds  to  evaporation  of  the  major  portion  of  the 
fuel,  is  always  found  in  a  fixed  region  of  the  chamber,  independent  of 
the  absolute  length  of  the  chamber.  The  experimental  results  appear  to 
be  In  agreement  with  the  theory  for  the  conditions  investigated. 


1. 


INTRODUCTION 

The  processes  involved  In  the  combustion  phenomena  in  a  liquid 
propellant  rocket  consist  of  spraying  of  fuel  and  oxidizer  into  the 
chamber,  mixing,  evaporation  of  the  dropletn,  inteidiffusion  of  the 
vapors,  and  then  the  chemical  reaction.  The  chain  of  events  stated 
here  is  still  an  oversimplified  version  of  the  true  account.  The 
overall  process  is  so  complex  that  ve  must  resort  to  the  study  of 
simplified  models  having  sane  of  the  most  important  features  of  the 
real  physical  phenomena.  In  the  evaporation  and  combustion  of  a 
spray,  we  simplify  to  a  point  where  vre  discuss  the  ballistics,  evapo¬ 
ration  and  combustion  of  individual  droplets  of  fuel.  If  ve  assume 
that  interaction  between  single  droplets  in  negligible,  then  evapo¬ 
ration  rates  of  a  distribution  of  droplets  are  additive.  Thus,  a 
distribution  of  droplets  can  be  used  to  approximate  the  evaporation 
rates  of  certain  spray  configurations. 

This  paper  depicts  combustion  in  a  rocket  motor  as  controlled  by 
the  evaporation  process  of  an  arbitrary  number  of  droplets  of  injected 
fuel.  The  evaporation  rate  of  the  droplet  of  fuel  is  of  course  most 
affected  by  its  immediate  environment,  which  in  this  case  happens  to 
be  high  velocity  gases  at  temperatures  in  the  range  of  5800°R,  and  at 
a  pressure  of  approximately  11  atmospheres.  The  environment  is  de¬ 
termined,  to  a  large  extent,  by  the  rate  of  evaporation  of  the  droplet 
of  fuel,  that  is,  by  the  rate  at  which  chemical  energy  is  added  to  the 
system.  Thus,  there  is  a  strongly  coupled  system  consisting  of  the 
evaporation  processes  and  the  fluid  dynamics  of  the  environment.  The 
results  of  the  evaporation  of  the  droplets  of  fuel  are  weighted  to  yield 
a  mass  rate  of  flow  which  is  consistent  with  approximately  a  500  pound 


thrust  chamber. 


a. 


THE  DROPLET  KOi.r-rL 

Tbo  evaporation  of  a  spray  is  treated  as  the  summation  of  the 
evaporation  of  individual  drops.  The  general  reliability  of  droplet 
drag  and  vaporization  rate  correlations  for  single  droplets,  when 
applied  to  sprays,  in  suspect.  This  work  is  aimed  at  determining  the 
consequences  of  the  assumption  of  additivity  of  droplets. 

— ",M-  'rh0  ]  ltlula  is  completely  mixed;  has  a  uniform  temperature, 
and  ite  composition,  vapor  pressure,  surface  tension,  specific  heat, 
and  viscosity  is  that  of  JP-5;  the  vapor  pressure,  surface  tension, 
specific  heat  and  viscosity  are  taken  a3  functions  of  temperature. 

— — ~  11116  fil3n  connate  of  two  components:  (a)  with  properciea 
of  JP-5A,  and  (b)  with  properties  averaged  between  carbon  dioxide  and 
water  vapor,  including  average  critical  pressure  and  temperature. 

Outer  Boundary  of  Film:  The  film  boundary  temperature  is  that  of 
the  bulk  gas;  where  the  temperature  and  pressure  of  bulk  gas  vary  in  a 
prescribed  manner  (from  the  hydrodynamic  calculations):  the  partial 
pressure  of  the  hydrocarbons  equals  zero  and  the  partial  pressure  of  the 
inert  component  equals  the  total  pressure. 

Bulk  Gas:  The  bulk  gas  temperature,  pressure  and  velocity  are  deter¬ 
mined  by  fluid  dynamic  calculations.  Its  composition  is  determined  by 
assuming  that  all  of  the  oxidizer  is  vaporized  and  that  chemical  equili¬ 
brium  exists  between  the  vr.porlzed  hydrocarbons  and  the  oxidizer. 

Inner  Boundary-  of  the  Film:  The  inner  boundary  temperature  is  that 
of  the  liquid  where  the  partial  pressure  of  the  hydrocarbon  equals  the 
vapor  pressure  and  the  partial  pressure  of  the  combustion  products  equals 
the  total  pressure  minus  the  partial 


pressure  of  the  hydrocarbon. 
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5£SL Trans ferj  The  heat  transfer  process  is  assumed  to  be  by  convection 
only  where  the  correlation  of  Ranz  and  Marshall  (l)  5  ig  uged> 


(Nu)h  r  ^  r.  1  +  0.6 


(1) 


r,  droplet  radius,  ft. 

h,  film  coefficient,  BTU,/ft2~sec-deg.  F 

Km,  mean  film  conductivity,  BTU/ft-sec-deg  F 

C  |i 

Pr,  Prandtl  number,  dimenaionless  =  -£ — 

K 

He,  Reynolds  number,  dimensionless  =  -^r  ( u**v ) 

d 

A  correction  is  made  for  sensible  heat  carried  by  the  evaporating 
fuel  moving  away  from  the  drop,  (6) 


hALCta-TL)2 


(2) 

Here  q^.  is  the  heat  in  BT'J  per  second  absorbed  by  the  droplet  by 

convection;  A_  is  the  area  of  the  droplet  in  square  feet,  and  T  ar^  T 

g  L 

are  the  gas  and  droplet  temperatui-es,  respectively.  Z  is  defined  by 


z  - 


p 


vy__c 

h  Al 


It  1b  Been  then,  that  Z  is  that  fraction  of  total  heat  transfer 
frern  the  gas  that  arrives  at  the  surface  of  the  liquid  drop.  (3  is  di¬ 
mensionless  since  W  is  the  vapor  from  the  evaporating  drop  in  pounds  per 
second,  and  is  the  specific  heat  at  constant  pressure  of  the  vapor  in 
KFU/lb  °F. 

Numbers  in  parentheses  indicate  references  at  end  of  paper. 


the  t-perrtuM  cbarEe  Bay  than  ha  cotputod,  slnon  x,  the  lats„t  hMt 
of  the  drop  (BTU/Ib)  to  to*.  If  e  to  the  time  In  .econde,  then 

.  nrnn 

I  ^  ®  CLML 


(3) 

where  qL  ,  the  net  heat  absorbed  by  the  evaporating  droplet,  iB  given 
by 

\  =  %  -  w  \ 

CL  is  the  specific  heat  of  the  liquid  (KCU/lb°F)  and  ^  ia  the  weight 
of  the  liquid  drop  in  pounds. 

Mass  Transfer  is  assumed  to  be  by  film  diffusion.  Here  thenaal 
diffusion  (the  Dufour  Effect)  is  neglected  and  the  correlation  of  Renz 
and  Marshall  (l)  is  used: 


w,=  =  2  *0.6  Sc’"  Re* 

(4) 

where  IJu  ,  Nusoelt  number,  dimensionless 
Dv  *  diffusivity,  ft2/eec 

Mg  »  molecular  weight  of  fuel  vapor,  lb/lb-mole 

R  ,  gas  constant,  ft-lb/lb-mole  °R 

Sc  ,  Schmidt  number,  dimensionless  -  — ti _ 
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The  dlffuaivlty  ia  calculated  by  the  equation  given  by  Bird  (2) 

■\  -  fcq  Tj'“  £  „ 

fl  Ma  KhS'i  p  (R 

VMa+  Mbj 


Here 


(5) 


molecular  diffusion 
coefficient 


Pca  iB  the  Pseudo-critical  pressure  of  JP. ■%,  atmospheres 
Pcb  ±a  the  average  critical  pressure  of  C02  and 

T  Z£?lr  t“PC”tUre  0f  tooplrt  «1  <X«ta*Mri 

Tca  is  the  Pseud-o-crltical  temperature  of  JP-5A,  °R 
Tcb  18  811  avera6e  critical  temperature  of  C02  and  H^O,  °B 
while  a'  la  a  correlation  constant  =  3.862  x  10^ 

004  ^  iB  a  correlation  constant  -  1.8229 

Ma  13  the  “oscular  weight  of  JP-5A,  lb/lb-mole 

"b  la  the  avera«e  molecular  weight  of  COg  and  H,,0,  lb/lb-mole 

A  correction  is  made  for  the  effect  of  bulk  motion  of  the  fluid  away 
fran  the  evaporating  interface  (3).  If  ^  le  the  vapor  pressure  of  the 
liquid  fuel,  and  p  the  total  static  pressure  of  the  gases,  then 


W  =  K3  AlPlck 


(6) 


6. 


where 


n 


_P 

P-pT- 


(7) 

a  is  a  correction  factor  to  account  for  the  unidirectional  diff¬ 
usion  process  rather  than  just  considering  the  equimolal  diffusion  trans- 
fer  coefficient  Kg  which  is  computed  fran  equation  k. 

The  change  in  radius  of  the  fuel  drop  in  time  can  be  calculated  if 
the  liquid  density  is  known  as  a  function  of  the  drop  temperature. 

The  continuity  equation  for  the  liquid  drop  is: 


d ( pV ) 
do 


=  -w 


v1  is  the  drop  volume 


•TT  r3 


4r-  W  r  f d  A  \  d  \ 
de  AL/>L  -3/>L[d XI  d0 


(8) 

Di^S  is  calculated  from  Ingebo's  correlation  on  liquid  drops  (k). 
The  drop  acceleration  is  known  once  the  gaB  density  J)  and  velocity  U 
is  stated. 


dV  _  3  r 

de  ~  "  a  0 


Z  (Lkvf 
A  ~  r 


(9) 

and  is  taken  to  have  the  same  sign  as  (U  -  V).  V  is  the  drop  velocity 
in  ft/sec  and  CQ  }the  drag  coefficient  is  given  by: 
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(10) 

Drop  breakup  was  considered  necessary  since  high  pressures,  above 
the  critical  temperature  of  the  fluid,  would  be  used.  Liquid  tempera- 
turea  could  not,  physically,  go  to  the  critical  without  induction  of 
aurface  tension,  <T  L  ,  to  the  point  where  droplet  stability  is  impaired. 
Therefore,  provision  was  made  for  the  droplet  to  be  replaced  by  an  arbi¬ 
trary  number  of  smaller  di-ops  of  equal  mass  whenever  on  arbitrary  Weber 
number  was  exceeded. 


We  rn  Lrju-v)\p 

<£ 

(11) 

Tliis  number  is  a  ratio  of  the  distorting  force  to  the  restoring  force. 

This  is  Intended  as  a  zero  order  approximation  to  qualitatively  detennine 
wliEvfc  efiecto  muy  occur. 


THE  HYDRODYNAMIC  MODEL 

The  thermodynamic  properties  of  the  combusted  fuel  and  oxidizer  were 
obtained  from  the  Second  Law  restriction  that  the  Free  Energy  must  be  a 
minimum  at  a  given  value  of  pressure  and  temperature  for  an  equilibrium 
condition.  The  reaction  equations  used  are: 

co2-  co  —  o  { co2  +in&o-^co~-  oh 

ICOi-ICO  —  Ol  CO  +  HaO  -  COi  — hz 

2  CO+f  M2O-|C02  — -  h 


0. 


These  equations  are  written  in  the  form  of 

AGj  =  2  A 

L  J 

(A) 

where 


and  d  ij  ls  the  coefficleat  of  the  i  th  component  of  the  j  th  reaction. 

is  the  number  of  moles  of  specie  i .  The  change  in  n^  species  are 
limited  "by  the  restraint 


(B) 

Here  £  io  the  extent  of  the  reaction. 

Also  at  a  fixed  T  and  p 


dG  =  X ML  dnt 

l 


A  parameter  X  is  introduced  such  that 


(C) 


=  -  A  Gj  d  A. 


Then  with  (A),  (b)  and  (c)  we  obtain 

dG  =  -[Z^G/|c|A=  lAdnL 

J  }  l 


(D) 
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If  d\>0,  then  a  thermodynamically  possible  change  in  G  can  occur,  i.e., 

G  decreases.  From  this  form,  the  sum  of  the  squares  of may  be  computed 
and  the  computation  terminated  when  ££  ;  £  is  determined,  frm 

the  accuracy  of  the  free  energy  data.  A  more  detailed  discussion  Is  given 
in  (?)• 

Once  the  minimum  value  of  free  energy  is  determined,  then  the  species 
n±  are  taaovn  and,  with  a  knowledge  of  the  specific  heat  data,  the  enthalpy 
and  molecular  weight  of  the  mixture  may  be  determined.  These  two  values 
are  used  in  the  energy  equation  (10)  and  the  equation  of  state  (12). 

The  relations  for  the  steady  state  equations  may  be  obtained  frcn 
the  nonsteady  equations  (5).  These  equations  relating  the  gas  pressure  P, 
density  p  ,  temperature  T,  and  velocity  U  are  the  equation  of  state: 

p  =  (tH/>  T 


_  (12) 

)  is  the  gas  constant. 

The  continuity  equation: 

„  d(fiUA)  =  w 
be  b  x 


(13) 

o 

A  is  the  chamber  cross  sectional  area,  ft 
W  is  the  mass  added  per  unit  length  per  unit  time. 

The  momentum  equation: 

iu  +  uau-.i  3p 
do  dx  p  dx  1 


(It) 


f  la  the  force  between  the  liquid  drops  find  gas, 


ft-lbf 

lbm 
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The  energy  equation  following  a  gas  particle  Is: 


d& 


U 


ah  _ 


<3>X 


c>p  ,  u  ap 


—  Cj  +  H 

1  p  be  p  ax 


(15) 

This  equation  can  be  reduced  to  a  relation  between  thermodynamic 
states  once  the  relation  between  enthalpy  (h)  and  temperature  is  known. 

For  the  steady  state  all  time  derivatives  may  be  dropped  from  (13) 
and  (l4)  so  that: 

ciCjdUA)  =  W  d* 

UdU=-  ~p  d  P  fd  X 


The  integrated  continuity  equation  becomes 


(p>JA\  = ;J~\ Wdx  =  p* 


The  integrated  momentum  equation  becomes 


(16) 


X+  AX 


—  P-*  + 


A 


Xi-/sy 


J 


U*  -  Ux*&x 


(17) 

The  energy  equation  for  a  control  volume  is: 

'  (h,„-  4>,)(CPTf  -p)  -  1 „ 


(18) 
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Here  q^  is  obtained  from  equation  (2). 

V  is  the  drop  velocity  and  Tf  1b  th-  flame  temperature. 


Equations  (12),  (l6),  (l7),  and  (l6)  CQffibined  ^ 
vritten  in  finite  difference  form,  from  the  set  of 


(3),  (6)  and  (9), 
equativns  to  be  solved. 


me  strategy  of  Computation 

The  calculation  was  divided  into  two  main  parts.  In  the  first  part 
a  single  droplet  is  followed  as  it  moves,  evaporates,  and  breaks  up,  in 
a  prescribed  environment  of  gas  pressure,  temperature,  and  velocity. 


Equations  (3),  (8),  and  (9)  are  simultaneous  ordinary  differential 
equations  in  the  variables  TL  ,  r,  and  v,  and  can  be  solved  if  initial 
conditions  plus  values  of  T,  p,  and  U  as  all  points  are  known.  This  is 

repeated  for  each  drop  size  present,  and  the  results  weighted  to  give 
the  total  evaporation. 


In  the  second  part,  the  hydrodynamic  equations  are  solved  for  a 
specified  evaporation  profile  taken  from  the  results  of  part  one. 
Equations  (12),  (l6),  (17)  and  (18)  are  solved  simultaneously,  since 
0  (x)>  \(x)>  f(x)  are  known  from  the  evaporation  calculation.  The 
boundary  conditions  for  the  gas  dynamic  equations  are  given  as: 


x  "  x0  P  =  P0>  T  =  Tf,  U  =  0 

at  x  =  L  UT  =  /2>ht 

i>  L 

100  basic  increments  were  subdivided  into  a  number  of  sub-steps 
which  were  determined  by  the  slope  of  the  evasion  curve.  With  10 
Bub-steps,  satisfactory  results  were  usually  obtained-  Normally,  1000 
Increments  were  used,  unless  the  liquid  temperature  changed  by  more  than 
50°F  in  an  increment,  in  this  case,  the  increment  size  was  reduced  by 


12. 


a  factor  of  10  for  that  sub-step.  If  the  liquid  temperature  chouse  vao 
still  too  large,  then  then  the  sub-step  was  reduced  by  102,  etc.  This 
insured  slowly  changing  integrands  and  provided  satisfactory  convergence 
for  all  cases  tried.  Physical  properties,  uhicta  were  needed  in  compu¬ 
tation  of  equations  (l)  thru  (11)  were  provided  by  a  series  of  subroutines, 
using  a  variety  of  equations,  as  well  as  tables  with  interpolations.  Change 
of  any  basic  equations  would  require  altering  the  evaporation  subroutine. 
Change  in  physical  properties,  other  than  input  quantities,  can  be  made 
by  replacing  the  appropriate  subroutine. 

In  the  evaporation  subroutine,  droplet  breakup  takes  place  according 
to  the  Weber  number  criterion,  equation  (ll). 

After  computing  evaporation,  the  main  program  weights  the  results  to 
compute  total  evaporation  at  each  Increment,  which  is  part  of  the  input 
information  to  the  fluid  dynamics  subroutine. 

In  equations  (l2),  (l6),  (17)  and  (l8)  form  the  heart  of  the  fluid 
dynamic  subroutine.  The  program  uses  information  given  to  it  by  the  e» 
vaporation  subroutine  to  produce  a  matrix  (3  x  100)  of  gas  pressure, 
temperature,  and  velocity.  That  is,  the  eonbustion  chamber  and  nozzle 
are  divided  into  100  sections  (l,  2,  3, . n,  n  +  1 . 100).  A 

separate  calculation  is  carried  out  in  each  section,  that  is,  going  from 
(n)  to  (n  +  1). 

The  area  which  is  available  to  gas  flow  A  ,  A  ,  Is  calculated 

n  n  +  1 

by  a  subroutine  program  and  also  the  total  mass  that,  is  present  at  any 
section  $5^  +  1  is  obtained  by  a  function  program.  The  subroutine 

checks  the  ratio  of  the  terns,  i-e-,  and  if  a  prescribed  limit  iQ 

exceeded,  the  step  size 


X 


nM 


">n 


n  A  X, 
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becomes 


ZlV  - 


NJ 


where  N  may  be  any  convenient  number  such  that 


h 


<Pn  *  1’ 

Here  the  primes  denote  the  new  coordinate.  The  ccmiputational  procedure 

is  carried  out  N  times,  then  the  ratio  is  checked,  and  so  on.  The 

method  of  solution  consists  of  assuming  a  T  ,  and  solving  for  U 

n  +  1  .  n  +  1 

in  equation  (l6).  Then  one  may  solve  for  p  knowing  the  U  from 

n  +  i  n  +  1 

equation  (17)  and  pn  from  ( 12 )  -  Tn  +  x  may  be  solved  for  in  equation 
(l8)  and  if 


T  -  T 

04i  1  n  1 1 


<  £ 


e  >o 


then  we  Bay  that  the  system  of  difference  equations  (12),  (l6),  (17),  and 
(18)  are  satisfied.  We  now  proceed  in  the  space  calculations  by  assuming 
TQ  +  2  and  continue  in  a  similar  manner-  The  superscripts  denote  iteration 
number.  Transfer  after  completion  of  100  steps  is  to  the  main  program, 
which  controls  the  drop  and  gas  calculation ■ 


DESCRIPTION  OF  EXPERIMENTAL  APPARATUS 

The  first  series  of  tests  were  conducted  with  liquid  oxygen  and 
JP-5A  in  a  two  inch  diameter  chamber,  employing  a  shoverhead  type  in¬ 
jector  and  a  simple  converging -diverging  nozzle,  (Fig.l).  The  only 
variable  in  this  first  series  of  tests  was  the  chamber  length,  i.e., 
using  the  same  in.ieotor-no?,zle  configuration,  and  the  same  propellant 
flow  rates,  the  chamber  length  was  to  be  varied  frcm  8  to  24  inches. 
Therefore,  the  chamber  may  be  composed  of  one  or  more  sections  depending 
upon  the  total  length  desired.  The  three  basic  components,  (injector, 
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chamber,  and  nozzle  ),  are  beld  to5ether  by  near,,  of  .  mechanical 
olanplns  am^ement.  AU  component.  a*,  uncoolea,  and  are  fabricated 
fro.  commercial  bra.e,  wth  the  exception  of  the  nozzle,  ttl,  l0  Bom. 
poeed  Of  a  graphite  liner  la  a  bran  collar,  and  la  alao  oncooled.  A 

complete  deecrlptlon  of  the  feed  eyztem,  co«t„l  eyetem,  end  rocket  motor 
can  be  found  in  (8). 

The  design  of  the  nozzle  for  the  steady  state  program  was  influenced 
by  the  future  objectives  of  the  Laboratory.  Among  these  objectives  was 
an  investigation  of  the  non  linear  aspects  of  combustion  instability.  In 
the  theoretical  analysis  of  combustion  instability,  the  effect  of  the  Mach 
number  of  the  gases  Is  highly  significant,  as  indicated  in  (3).  Briefly 
stated,  if  (l-M2)  is  approximately  unity,  i.e.,  if  K  ie  approximately 
0  2  or  less,  the  mathematical  analysis  can  be  greatly  simplified  by  means 
of  linearizing  assmptlons.  However,  the  actual  problem  is  non  linear, 
and  in  practice  the  Mach  numbers  at  the  nozzle  entrance  are  well  above  the 
"linear"  value  of  0.2.  Since  the  entire  steady  state  program  was  being 
conducted  as  a  preliminary  step  to  a  non-steady  non  linear  analysis  of 
combustion  instability,  the  nozzle  contraction  ratio  was  designed  for  a 

nozzle  entrance  Mach  number  of  .1*5,  based  upon  isentropic  flow  and  constant 
ioentropic  exponent. 

The  verification  of  the  nozzle  entrance  Mach  number  required  a 
measurement  of  the  total  pressure  at  that  point.  This  entails  designing 
a  probe  that  can  exist  in  an  extremely  high  temperature,  oxidizing  atmos¬ 
phere.  A  water-film  cooled  graphite  probe  was  developed  and  employed 
successfully,  (Fig. 2).  The  cooling  water  ij  sprayed  through  the  probe 
body  into  the  gas  stream  at  a  pressure  slightly  in  excess  of  the  gas 
pressure.  Thus,  a  film  type  cooling  is  obtained.  The  mass  flow  of  water 
Into  the  chamber  amounts  to  only  1  %  of  the  propellant  flow. 
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The  Instrumentation  necessary  for  the  measurement  of  Btatlc  pressure 
gradients  in  a  rocket  motor  was  selected  after  a  consideration  of  several 
factors:  namely;  proximity  of  measurements,  expected  pressure  gradient, 
accuracy,  and  cost-  Assuming  tnat  the  pressure  was  to  "be  measured  at 
eight  longitudinal  locations,  the  obvious  solution  would  be  to  install 
eight  transducers  down  the  length  of  the  chamber.  However,  each  of  the 
considerations  above  eliminated  this  technique.  Tests  conducted  on  the 
shorter  length  chambers  would  require  a  pressure  measurement  at  one  inch 
intervals.  The  physical  size  of  the  diaphragm  on  a  water  cooled  trans¬ 
ducer  varies  from  l/2  to  1  inch  depending  on  the  manufacturer.  Therefore, 
each  transducer  would  actually  measure  an  average  pressure  over  an  inter¬ 
val  comparable  to  the  center  to  center  distance  between  pressure  taps. 

In  addition,  it  was  expected  that  the  gradient  at  certain  sections  of  the 
motor  would  be  as  small  as  1  psi  in  2C0  psi  between  adjacent  taps.  UBing 
a  500  psi  full  range  transducer,  a  drift  of  .1  °f  of  full  scale  in  oppo¬ 
site  directions  by  two  adjacent  transducers  would  yield  zero  pressure 
gradient.  Rather  than  develop  a  method  of  compensating  for  the  different 
drift  characteristics  of  adjacent  transducers  and  associated  electronics, 
it  appeared  more  feasible  to  develop  a  technique  that  employed  a  single 
transducer  to  measure  the  pressure  gradient.  In  addition,  if  such  a 
system  could  be  developed  the  dollar  savings  to  the  experimental  program 
would  be  significant.  Rather  than  employ  eight  water  cooled  transducers 
it  would  be  possible  to  use  a  single  transducer.  It  will  subsequently  be 
shown  that  the  transducer  need  not  be  cooled,  restating  in  an  additional 
savings.  All  of  the  above-mentioned  factors  lead  to  the  development  of 
a  pressure  scanner  (Fig. 3)- 
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The  scanner  permit.  .  single  transducer  to  travel  -fra  pressure  tap 
to  pressure  tap  and  consecutively  record  tie  pressure  at  eight  different 
location,  in  the  chanter.  The  pressure  scanner  consist,  of  three  piston 
and  cylinder  units.  Units  A  and  B  each  contain  facilities  for  measuring 
the  pressure  at  eight  locations,  while  unit  c  is  a  hydraulic  driver.  The 
taps  located  in  the  rochet  rotor  are  extended  to  the  bulkhead  of  either 
unit  A  or  B  and  then  to  compartments  formed  by  the  Inner  oylinder  or 
piston,  V  rings,  and  the  outer  travelling  cyiinder  which  house,  a  trans¬ 
ducer.  The  reciprocating  motion  of  the  cuter  cylinders,  (which  enables 
the  transducer  to  consecutively  measure  the  pressure  at  adjacent  tap,  in 
the  rocket  motor)  is  governed  by  a  double  acting  hydraulic  piston  and 
cylinder;  unit  0.  Previsions  have  been  made  for  varying  the  speed  and 
length  of  stroke  of  the  unit,  In  addition,  a  method  of  continuously 
recording  the  erect  location  of  the  transducer  housing  has  bee»  developed. 
This  consists  of  measuring  the  voltage  drop  between  ,  moving  electrical 
contact  mounted  on  the  transducer  housing  and  a  fiat  copper  bar  mounted 
alongside  the  scanner  cylinder.  The  copper  bar  is  actually  an  eight  step¬ 
ped  voltage  divider.  The  length  of  each  step  is  equal  to  the  length  of 
the  cmnpartments  in  the  scanner  cylinder.  The  distnnee  between  steps  or 
sections  is  equal  to  the  thickness  of  the  "0"  rings  separating  the  compart¬ 
ments.  Each  of  the  copper  bar  sections  is  a  a  different  predetermined 
ag  Therefore,  by  simultaneously  recording  the  output  of  the  trans¬ 


ducer  and  the  voltage  drop  to  the  moving  electrical  contact,  it.  is  possible 
to  match  the  pressure  magnitude  with  the  location  at  which  it  was  measured. 

The  output  of  the  transducers  is  recorded  on  a  recording  oscillograph. 

To  record  the  absolute  pressure  at  each  port  would  result  in  reduced  resolu- 

Therefore,  only  the  pressure  difference  about  a  fixed  point  is  recorded. 
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Thus  kQ  psi  rather  than  200  Pst  is  full-scale  deflection  on  the  recorder. 

The  above  technique  is  accomplished  by  several  methods  depending  upon  the 
type  of  transducer  employed  in  the  scanner.  Using  a  piezoelectric,  crystal 
transducer,  the  crystal  is  grounded  until  the  starting  transient  has  been 
completed .  The  pressure  on  the  transducer  diaphragm  at  the  instant  that 
the  ground  connection  is  broken  and  all  subsequent  pressures  to  be  measured 
is  recorded.  To  determine  the  absolute  pressure  at  each  port,  from  the 
records  of  the  gradients,  requires  a  measurement  of  the  absolute  pressure 
at  one  point  in  the  chamber. 

Several  objections  may  arise  as  to  the  accuracy  and  applicability  of 
the  above  described  system  for  the  measurement  of  pressure  gradients  in  a 
rocket  motor.  It  is  obvious  that  the  pressure  at  each  chamber  tap  is  not 
measured  simultaneously,  but  rather  there  is  a  time  lag  between  the  measure¬ 
ment  at  successive  ports.  However,  the  speed  of  the  scanner  can  be  adjusted 
so  that  a  complete  traverse  of  the  chamber  can  be  made  in  one  second.  Since 
a  normal  tent  firing  lasts  Bix  seconds  after  steady  state  conditions  have 
been  reached,  six  complete  traverses  of  the  motor  can  be  obtained  in  one 
run.  Examination  of  the  data  indicates  that  one  traverse  is  identical  to 
the  next,  both  as  to  gradient  and  absolute  pres cure  at  each  port.  Thus,  it 
1b  concluded  that  once  the  maximum  chamber  pressure  has  been  reached,  the 
motor  operation  is  steady  in  time  unless  an  instability  arises. 

In  spite  of  the  previous  remarks  regarding  tho  steadiness  of  operation 
and  the  reproauclbltlty  of  gradients  on  a  recording  trace  during  a  te3t  run, 
the  pressure  at  each  chamber  location  is  not  invariant.  Under  the  smoothest 
conditions  there  is  Btiil  a  good  deal  of  combustion  noise  present.  In  ad¬ 
dition,  as  pointed  out  in  (3)  and  (U),  there  exists  an  instability  regime 
for  any  particular  oxidizer-  fuel  ratio  of  a  given  set  of  propellantB. 

During  the  course  of  the  tests,  "rough"  or  "unstable"  operation  wsb  encountered 
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on  several  occasions,  particularly  with  the  shorter  chanters,  ranging 
in  length  from  eight  to  twelve  inches.  Further  discussion  of  the  oc¬ 
curence  of  pressure  oscillations  in  excess  of  combustion  noise  is  pre- 
sented  in  the  discussion  of  experimental  results.  Despit  _*  the  above 
comments,  it  is  still  possible  to  determine  the  steady  state  pressure 
gradients,  where  steady  state  is  defined  as  the  average  or  D.C.  level 
of  the  transducer  output. 


EXPERIMENTAL  RESULTS  -  COMPARISON  WITH  THEORY 

Steady  state  pressure  gradients  were  measured  in  a  two  inch  diameter 
variable  length  rocket  motor.  Initially,  static  pressure  was  measured  at 
seven  longitudinal  locations  and  total  pressure  was  measured  approximately 
one  inch  from  the  nozzle  entrance.  After  it  was  determined  that  the  nozzle 
entrance  Mach  number  was  greater  than  0.4  the  total  pressure  measurements 
were  abandoned.  Subsequently,  it  became  possible  to  measure  the  static 
pressure  gradient  at  eight  longitudinal  locations. 

The  first  series  of  tests  were  conducted  on  a  17  inch  chamber  length - 

The  results  are  plotted  in  Figure  (4)  Results  indicate  the  prescenee  of 

a  steep  pressure  gradient  between  four  and  ten  inches  from  the  injector 

face.  The  data  indicates  that  a  29  psl  gradient  was  measured  between  1  i 

3 

and  14  g  inches  from  the  injector  face.  This  amounts  to  17%  of  the  pressure 
measured  at  the  1  ^  inch  tap.  In  the  region  between  four  and  ten  Inches 
from  the  injector,  (35%  of  the  chamber  length),  a  gradient  of  21  psi  was 
measured.  This  is  equal  to  72  %  of  the  measured  total  gradient,  or  a  some¬ 
what  smaller  percentage  of  the  gradient  obtained  by  extrapolating  to  the 
injector  face  and  nozzle  entrance.  According  to  the  analysis  set  forth  in 
(l),  the  major  portion  of  combustion  takes  place  in  this  steep  pressure 
gradient  region.  The  gradient  in  the  remaining  portion  of  the  chamber. 
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(10  -  17  inches),  is  flat;  7  psi  in  7  inches.  It  is  felt  that  this 
portion  of  the  pressure  drop  is  due  to  heat  transfer  nnrj  low  friction. 

The  ratio  of  static  to  total  pressure  is  0-915  at  the  It  g  inch  tap.  This 
corresponds  to  a  Mach  number  of  0.41  for  a  gas  whose  laentroplc  exponent 
is  1.11.  The  value  of  1.16  waa  determined  by  the  analysis  set  forth  ini 
(l).  The  oxidizer  fuel  ratio  for  this,  and  all  subsequent  tests  discussed 
in  thiB  report  is  2.76  +  5%  •  The  total  propellant  flow  ia  2.47  lbm/sec 
45^.  The  O/F  equivalence  ratio  is  0  -  '793  •  The  percentage  deviation 
Btated  above  does  not  account  for  the  experimental  error  in  a  single  test, 
but  1b  due  to  a  variation  in  measured  flow  from  one  teat  to  another. 

The  theory  predicted  that  the  main  evaporation  and  combustion  of  the 
drops  oecured  within  a  fractional  section  of  the  chamber  length.  Here  the 
fluid  dynamic  flow  field  increased  the  evaporation  rate  which  in  turn  in¬ 
creased  the  velocity  field  and  hence  the  pressure  gradient. 

Typical  results  of  the  integration  of  equations  (3),  (8),  (9),  (12), 
(l6),  (17)  and  (l8)  are  shown  in  Fig. 4.  Curve  A  is  the  solution  for  p(x) 
when  the  boundary  condition  for  temperature  at  the  injector  is  equal  to 
the  flame  temperature  (5877°^).  Curve  B  is  the  solution  for  p(x)  when  the 
gas  temperature  at  the  injector  face  1b  5500°R-  Curve  B  is  intended  to 
show  the  effect;  of  heat  transfer  or  combustion  efficiency  on  the  pressure 
drop.  Curve  B  allows  a  greater  masB  flow  through  the  nozzle  due  to  lower 
chamber  temperatures.  <pma5(  =  2.36  is  the  experimentally  measured  flow  rate. 

In  order  to  determine  the  effect  of  chamber  length  on  the  position  of 
the  steep  portion  of  the  gradient,  and  hence  the  combustion  process,  the 
chamber  was  reduced  to  8  g  inches  with  subsequent  tests  to  be  ruii  with 
chamber  length  increments  of  two  to  three  inches-  The  gradients  obtained 
in  an  8  g  inch  chamber  and  a  10  g  inch  chamber  are  presented  in  Figure  5. 
Both  curves  indicate  a  continuously  increasing  gradient  over  a  region  that 
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extends  approximately  one  inch  from  the  injector  face  to  one  inch  from 
the  nozzle  entrance ■  It  was  stated  previously  that  the  gradient  in  the 
IT  inch  chamber,  Figure^),  started  to  decrease  ten  inches  frem  the  in¬ 
jector  face.  As  a  result  of  the  Increasing  gradient  at  the  nozzle  en¬ 
trance,  it  is  concluded  that  evaporation  and  combustion  of  the  fuel  is 
not  completed  in  the  shorter  chambers-  This  is  in  agreement  with  the 
predictions  of  the  theoretical  analysis  which  is  also  shown  in  Figure  (5) 
and  Figure  (6).  The  Mach  number  at  a  point  one  inch  from  the  nozzle  en¬ 
trance  of  the  8  jj  inch  chamber  was  0.23,  as  compared  with  0.4l  for  the 
17  inch  chamber.  The  low  Mach  number  is  attributed  to  the  low  gas  gener¬ 
ation  rate  in  the  shorter  chambers.  In  addition,  the  velocity  increment 
of  the  gases  through  the  last  inch  of  the  chamber  is  much  greater  for  the 
shorter  lengths  than  for  the  17  inch  motor.  This  is  due  to  the  fact  that 
evaporation  and  combustion  is  still  present  in  the  short  chambers,  while 
it  virtually  has  been  completed  in  the  longer  chambers.  Results  fretn 
additional  tests  with  the  8  ^  inch  motor  Indicate  that  is  is  possible  to 
increase  the  Mach  number  by  increasing  the  O/F  ratio.  Since  the  main  body 
of  testB  were  run  with  fuel  rich  mixtures,  an  increase  in  the  O/F  ratio 
brings  it  closer  to  stoichiometric  proportions,  causing  a  more  complete 
combustion  process  in  the  chamber,  and  a  higher  Mach  number  at  a  point  one 
Inch  before  the  nozzle. 

The  experimental  results  obtained  with  a  12  J  ,  15  i  ,  and  19  inch 
chamber  are  reproduced  in  Figure  (6).  All  three  tests  indicate  a  decreasing 
gradient,  or  point  of  inflection,  approximately  ten  inches  from  the  injector 
face.  The  gradient  in  the  region  of  the  injector  face  appears  to  be  steeper 
in  the  case  of  the  12  ^  inch  chember,  than  in  the  longer  lengths.  This 
would  indicate  a  rapid  breakup  of  the  injected  spray,  and  hence  evaporation 
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and  combustion  closer  to  the  injector  face. 

The  above  phenomenon,  (steep  pressure  gradients  at  the  injector  face) 
is  also  present,  though  to  a  somewhat  lesser  extent,  in  the  8  ^  and  10  g 
inch  chambers.  The  steeper  gradients,  indicative  of  jet-breakup  and  cLrqp- 
i.ct  evaporation,  is  attributed  to  "rough"  or  "unstable"  operation.  Due  to 
the  low  frequency  response  of  the  recording  oscillograph  it  is  not  possible 
to  quantitatively  describe  the  pressure  oscillations.  However,  based  upon 
the  relative  amplitude  of  the  pressure  oscillations  and  the  frequency  of 
chamber  and  nozzle  burnout,  the  12  $  inch  chamber  would  be  the  most  un¬ 
stable  configuration.  The  oscillations  in  the  shorter  chambers  were  of  a 
lesser  degree,  whereas  the  longer  lengths  exhibited  only  normal  combustion 
noise.  It  is  felt  that  the  increased  pressure  oscillations  in  the  shorter 
chambers  causes  a  more  rapid  breakup  of  the  injected  spray,  resulting  In 
a  shift  of  the  combustion  zone  towards  the  injector  face.  Previous  work 
dune  at  the  Propulsion  Research  Laboratory,  (ll)  discuaseB  in  detail,  the 
effect  of  chamber  length  and  pressure  oscillations  on  liquid  jet  breakup. 

It  might  be  mentioned  at  this  time  that  the  steady  state  analytical 
work  did  not  include  the  .jet-breakup  distance.  It  appears  to  be  a  critical 
factor  for  the  initiation  point  of  the  evaporation  computation.  A  fixed 
of  2  inches  was  used  Id  all  computations  for  qualitative  purposes - 
A  more  complete  analysis  would  incorporate  a  mechanism  for  jet  breakup 
and  drop  formation.  It  would  appear  from  the  results  obtained  under  stable 
conditions,  that  this  2  inch  distance  is  of  the  right  order  of  magnitude. 
Iadiead  of  Rocketdyne  in  private  caramunication  has  also  indicated  a  2  inch 
region  with  similar  injector  configurations.  He  obtained  this  result  with 
streak  photography . 


A  question  may  arise  as  to  a  reason  for  the  variation  of  absolute 
uressure  at  the  injector  face  for  the  test  results  presented.  Actually, 
the  variation  is  only  3  per  cent,  and  is  equivalent  to  the  error  associated 
with  the  bourdon  tube  pressure  gauge  used  to  measure  the  reference  absolute 
pressure  at  a  point  in  the  chamber.  It  was  previously  stated  that  the 
reference  pressure:  was  used  only  to  convert  the  pressure  gradient  re¬ 
cordings  into  absolute  pressure  distributions.  The  significance  of  the 
variation  is  very  slight,  lu  should  be  stressed  that  the  main  function  of 
the  experimental  program  was  to  determine  the  pressure  gradient  for  a  given 
set  of  injector  boundary  conditions.  A  variation  of  1  or  2  psi  in  the  ab¬ 
solute  pressure  would  have  little  or  no  effect  on  the  gradient  for  a  given 
set  of  injection  parameters. 

Figure  T  shows  the  effect  of  initial  drop  radius  on  the  position  of 
the  pressure  gradient  in  the  15*5  inch  chamber.  The  three  drop  model  uses 
a  Rosin-Rommler  distribution  ( 12 )  to  describe  the  spray  distribution.  The 
small  drops,  which  make  up  ?8  °/Q  of  the  spray  (represented  by  the  26u  drop), 
evaporate  and  combust  quickly.  The  pressure  gradient  is  thus  initially 
steep  and  results  in  accelerating  the  evaporation  process  of  the  larger  drops* 
The  80  micron  drop  is,  thus,  comp] etely  evaporated  within  the  chamber.  On 
the  other  hand,  the  100  micron  model  evaporates  very  much  more  slowly  due  to 
the  absence  of  this  interaction,  and  so  a  15*5  inch  chamber  is  insufficient 
for  -proper  combustion. 

Figure  9  shows  the  gas  temperature  and  velocity  as  a  function  of  position 
in  the  15»5  inch  chamber*  Figure  10  shows  the  results  of  the  solution  of  the 
differential  equations  of  the  liquid  subsystems  for  the  1^*5  inch  chamber. 

The  liquid  temperature  rises  rapidly  but  then  remains  practically  constant 
throughout  evaporation*  ^/r  initially  Increases  due  to  expansion  and  then 


decreases  due  to  evaporation.  In  this  ca3e  r  ^  .19  x  lo”"*  ft,  V  =  1 

O  *”  '  c 

Ft/sec  and  Tq  =  600°R. 


CONCLUSIONS 

The  picture  is,  then,  one  of  slowly  evaporating  drops  flowing  down 
a  field  of  gradually  increasing  velocity  until  a  co-operative  effect  be¬ 
tween  the  drops  and  flow  field  causes  an  acceleration  of  evaporation.  In 
certain  cases  it  has  been  found  that  if  the  velocity  gradients  are  really 
severe,  shattering  of  the  droplets  occur.  This  happened  in  the  nozzle  of 
the  6  inch  chamber* 

When  the  Weber  number  is  given  a  value  of  100,  droplet  shattering 
(Fig.6)  occurs  in  the  chamber  (dashed  curve)  and  drastically  changes  the 
picture  of  evaporation.  Just  as  in  the  case  of  drops  in  a  spray  with 
continuous  evaporation  a  similar  effect  occurs  with  breakup.  When  the  drops 
start  breaking  up  evaporation  increases  very  quickly  causing  a  pressure  drop 
with  a  corresponding  increase  in  the  gas  velocity.  The  resulting  velocity 
gradient  then  increases  the  breakup,  etc.  The  basic  difference  between 
continuous  evaporation  model,  and  the  breakup  model  is  the  thickness  of  the 


For  a  given  set  of  injection  parameters  under  stable  operating  conditions, 
the  results  indicate  that  the  steepest  portion  of  the  gradient,  which  corre¬ 
sponds  to  evaporation  of  Die  major  portion  of  the  fuel,  is  always  found  in  a 
fixed  region  of  the  chamber,  independent  of  the  absolute  length  of  the  motor. 
For  the  tests  diBcussed  in  this  report,  the  steepest  portion  of  the  gradient 
was  located  between  four  and  ten  inches  :'rom  the  injector  face.  Chambers 
shorter  than  ten  inches  exhibited  combustion  in  the  nozzle.  It  was  also 
found  that  "rough"  or  "unstable"  operation  shifted  the  gradient.  Increasing 
pressure  oscillations  caused  a  shift  of  the  steep  portion  of  the  gradient 


towards  the  Injector  face* 

In  conclusion,  it  appears  that  evaporation  is  a  proper  rate  control¬ 
ling  mechanism  for  the  cases  investigated.  By  proper  shaping  of  the  combus¬ 
tion  chamber  the  evaporation  schedule  may  he  extended  ovci-  the  entire  chamber, 
length,  i.e.,  a  pear  shape  may  distribute  the  energy  release  pattern  (<J) 
evenly  over  the  entire  chamber.  An  inverse  effect  may  be  created  by  intro¬ 
ducing  constrictions  in  the  combustion  chsmber  to  localize  combustion  by 
shattering  the  droplets  in  a  high  velocity  gradient  field. 

Finally,  it  must  he  observed  that  the  intent  of  the  experimental  program 
was  to  determine  steady  state  pressure  gradients  for  verification  of  the 
theoretical  analysis  and  any  conclusions  based  upon  the  experimental  program 
alone  requires  further  substantiation.  The  analytical  work  will  be  extended 
to  a  nonlinear  nonsteady  instability  analysis  U3ing  the  results  of  the 
present  investigation. 
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